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HOTROEUCTIOW 

This  investigation  treats  a  truncated  sequential  test  for  testing 
a  sinple  hypothesis  against  a  siiitple  alternative.  Attention  is  con¬ 
fined  to  the  case  of  sacpling  froo  a  Poisson  population.  Although 
Wald's  /"ll?^  sequential  prohahlllty  ratio  test  (SPRT)  will  termi¬ 
nate  with  probability  one  at  some  stage  of  an  experiment,  there  is 
no  guaranteed  upper  bound  for  the  sample  size  of  this  test.  However, 
for  cases  in  which  time  and  cost  of  saE5>ling  are  involved  we  often 
have  to  face  the  situation  of  mal'wlng  a  definite  decision  within  a 
given  number  of  trials.  With  this  in  mind,  several  sequential 
tests,  such  as  Anderson's  Armitage's  /”27  and  Donnelly's  /”37» 

have  been  proposed  which  are  truncated  within  some  fixed  stage.  How¬ 
ever,  none  of  these  are  diirectly  applicable  to  Poisson  sampling. 

Here,  as  a  variation  of  the  SPRT,  we  treat  Hall's  minimura 
probability  ratio  test  (lIPBT).  ^Aiereas  for  the  SPRT  the  rejection 
and  acceptance  lines  are  parallel,  for  the  IIPRT  these  t«fo  lines  con¬ 
verge  so  that  the  test  always  terminates  with  a  definite  decision  by 
a  predetermined  stage. 

More  specifically,  we  present  and  coopare  several  test  procedures 
for  testing  idiether  a  Poisson  parameter  has  value  or  (speci¬ 
fied  numbers)  based  on  a  sequence  of  Independent  observations  from 
a  Poisson  population.  The  test  procedures  considered  here  are: 


The  numbers  in  square  brackets  refer  to  the  bibliography  listed  at 
the  end. 
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l)  the  sequential  probability  ratio  test,  both  untruncated 
(SPRT)  and  truncated  (SPRT^) (Chapter  I), 

ii)  the  niniimm  probability  ratio  test  (MPRT)  (Chapter  II ),  and 
iii)  the  most  powerful  fixed  saniple  size  test  (fsst)  (Chapter  IIl)> 

g 

A  diagram  for  carrying  out  these  tests  appears  in  Figure  I  (A. 2). 

A  requirement  in  all  of  these  test  procedures  is  that  each  error 
probability  should  not  exceed  a  common  specified  level  o.  The  bases 
for  comparison  of  the  test  procedures  are  the  operating  characteristic 
(OC)  function^  the  expected  sajiqple  size  (ASlf)  function^  and  the  stan¬ 
dard  deviation  of  the  sample  size  (SDR)  function.  Hoeffdlng's  lower 
bound  £’£!  (HUB)  on  the  ASH  at  an  intermediate  value  is  also  presented. 
(Chapter  III) 

Major  attention  is  given  to  the  MPBT  since  the  other  procedures 
are  quite  well  known.  An  extttislve  discussion  of  this  test  pro¬ 
cedure  appears  in  Chapter  II.  Actually,  the  MPRT  can  only  guarantee 
an  upper  hound  on  an  average  of  the  two  error  probabilities.  Achieving 
equal  error  probabilities  hinges  on  the  choice  of  an  intermediate 
value.  is  here  chosen  to  make  equel  the  divergence  (a  con¬ 
cept  from  information  theory)  between  and  and  between 
and  Xg  .  This  value  of  X^  is  denoted  D.  Some  consideration  is 
also  given  to  x^  •  s,  the  slope  of  the  SPRT  acceptance  lines.  The 
calculations  indicate  that  either  of  these  choices  is  quite  success¬ 
ful,  D  being  slightly  better  for  moderate  or  large  a  values  and 
8  being  sli^^tly  better  for  small  a  values  (<  .01).  (See  the 

tables  and  Chapter  V. ) 

2  — 

All  figures  and  tables  appear  in  the  Appendix. 
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For  the  IIPRT,  calculation  of  seven  points  on  the  OC  curve,  ASI^ 
curve  and  SDH  curve  were  carried  out,  a  »  .1,  .05,  .01  and  .001,  for 
the  following  pairs  of  hypotheses: 


Table 

I 

II 

III 

IV 

V 

VI 


.1  .5  D 

.5  .8  D 

.5  1.  D 

5.  8.  D 

.1  .5  B  e.nd  s 

1.  2.  D  and  s 


Also  presented  In  Table  I  and  VI  are  the  sari^ple  size  for  the 
FSST,  the  maxLiiiuu  sample  size  n^  of  the  MPBT,  HOeffdlng's  lower 
bound  at  D  and  s,  and  the  ASN  of  the  SPRT  at  four  values,  0, 

8,  and  Xg,  calculated  by  Wald's  approxloatloi  The  exact  cal- 

5 

culatlons  were  performed  on  the  DNIVAC  1105  by  a  program  described  In 
Chapter  IV  and  A.?. 

For  purposes  of  cooiparlng  the  three  test  procedures,  the  calcu¬ 
lation  of  l6  points  on  the  OC  curve,  ASH  curve,  and  SW  curve  were 
carried  out  for  all  tests,  for  a  »  .1,  .05,  .01,  and  .001,  and  for 
the  following  pairs  of  hypotheses: 

Table  Figure  ^1  ^2 

VII  2  .01  .1 

VIII  5  2.  4. 


5 


In  this  thesis  the  computer  refers  to 


UWIVAC  1105. 
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Exact  calculatkUB  irere  done  for  the  OC,  ASlf  and  SDN  fiuictlons  of  the 

IIPBT  (\  =  D)  and  the  SPRT  truncated  at  n^>  the  maximum  sample  size 
o  o 

of  the  I'lPRT,  and  Hoeffdlng's  lower  hound  at  \  ■  b  and  D  by 
the  computer.  Also  calculated  were  three  values  of  Niald's  approxima¬ 
tions  to  the  OC  and  ASN  functions  for  the  untruncated  SFBT.  The 
FSST  calculation  are  based  on  Poisson  tables  to  the  extent  avail¬ 
able  and  otherwise  on  the  normal  approximation.  A  discussion  of  the 
results  of  the  calculations  appears  In  Chapter  V. 

One  use  for  tests  concerning  a  Poisson  population  Is  In  quality 
c(xitrol  work  where  the  defects  In  a  unit  are  counted  and  the  quality 
of  a  lot  or  process  Is  Judged  by  the  average  number  of  defects  per 
unit.  This  differs  from  the  test  where  each  unit  Is  placed  Into 
a  "defective”  or  "non- defective"  category  and  the  quality  of  a  lot 
or  process  Is  determined  by  the  total  number  of  defectives. 

The  defects  per  unit  analysis  Is  useful  under  the  following 
conditions . 

I)  If  almost  every  unit  contains  at  least  one  defect,  a  fracti^ 
defective  plan  . . .  "defective"  or  "non-defective"  classification  . . . 
obviously  not  feasible  for  such  a  case. 

II)  For  products  which  are  e:qpenslve  to  produce  or  Inspect  and 
products  customarily  Inspected  in  small  lots,  It  is  too  costly  to  ob¬ 
tain  sai:q>les  large  enou^  to  assure  hif^  discrimination  by  a  fraction 
defective  plan.  However,  if  the  number  of  samples  observed  Is  suffi¬ 
ciently  large  and  the  quality  of  the  product  Is  high,  the  result  will 
not  differ  greatly  between  these  two  tests. 

III)  For  this  test  to  apply  exactly,  the  defects  must  be  randomly 
and  independently  distributed. 


CHAPTER  I 


SEQUENTIAL  PROBABILITY  RATIO  TEST 
AND  ITS  APPLICATION  TO  THE  POISSON  DISTRIBUTION 

1.1.  Ual&*s  Sequential  Probeblllty  Ratio  Test  (SPRT)  for 
Testing  a  Simple  Hypothesle  against  a  Simple  Alternative  /il7. 

Suppose  x^,  ....  is  a  sequence  of  Independent  observa^ 
tlons  with  a  cotmnon  density  function  f(x:d)  and  let 

x^j  ...  x^),  m  =  1,  2,  ...,  be  a  sample  of  size  m. 

Let  and  ot^  be,  respectively,  the  desired  probability  of  accept¬ 
ing  the  alternative  hypothesis  Hgtd  »  8^  when  the  true  parameter 
Is  6^,  and  of  accepting  the  null  hypothesis  H^:@  s  6^  when  the 
true  parameter  Is  Sg,  and  call  (a^,  citg)  the  strength  of  the  test. 
We  also  denote  by  dj^  the  decision  to  accept  Hj^(i  *  1.  2). 

For  any  positive  Integer  m,  the  Joint  density  function  f^^ 
of  a  sample  of  size  m  under  H^  (l  »  i.  2)  Is 

(1.1.1)  fjjjj  -  f(Xj;ep  .  f(x2:e^)  ...  f(x^:e^) 

Then,  the  SPRT  Is  carried  out  as  follows: 

For  suitably  chosen  A  and  B  (0  <  B  <  1  <  A),  at  the  m>th 
stage  (m  3  1.  2  ...)  of  the  experiment, 

(I)  stop  sampling  and  accept  H^  If 

(1.1.2)  <A  (j-1.2.  .hhI)  and  ® 

(II)  stop  sampling  and  accept  Hg  If 
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(1.1.5)  <A  (J=1.2...m.l)  and  fgm/^ln^  ^ 

(ill)  otherwise j  continue  sanqpllng  until 
either  category  (l)  or  (ll). 

The  calculation  of  A  and  B  to  obtain  the  desired  strength 
iot^,  ag)  Is  very  laborious.  Therefore,  In  practice  Wald  suggested 
putting 

(1.1.4)  A'5  (1  -  a^)/a^  and  B's  Og/d  -  a^^) 
as  substitutes  for  A  and  B,  respectively. 

Denoting  the  resulting  error  probabilities  by  and  we 
can  easily  see  that 

(1.1.5)  aj,+a^<a^+ag 

and  therefore,  at  least  one  of  the  Inequalities,  and 

ot^  <  ag,  must  be  satisfied.  Moreover,  If 
Inequalities  are  almost  achieved. 

The  most  Important  characteristics  of  the  SFBT  are  the  opera¬ 
ting  characteristic  (OC)  and  the  Aveirage  Sample  Nuniber  (ASN)  func¬ 
tions.  The  OC  function,  L(9),  Is  defined  as  the  probability  of 
accepting  the  null  hypothesis  when  the  true  parameter  Is  0. 

An  approximation  la  given  by 

.hC®)  1 

(1.1.6)  L(e)  -  ^h(e;  [ght®; 

where  h  satisfies 

(1.1.7)  Ee^=  1 

with  z  =  log  f(x:9g)/f(x:ej^)  .  If  h  -0  Is  the  only  solution  of 

(1.1.7) ,  the  right  hand  side  of  (1.1.6)  Is  evaluated  by  taking  the 
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limit  using  1' Hospital's  Rule. 

The  ASN  function,  £0(0)  >  1b  the  expectation  of  the  sample 
size  when  @  Is  the  true  parameter  and  Is  approximately  given  /Il_7 
hy 

(1.1.8)  ElgCn)-  |L(e)log  B+^L(e27log  Aj./Eg(z)  If  Eg(z)/ 0 
and 

(1.1.9)  Eg(n)-  |L(e)(log  B)®+2i-L(el7(log  A)2j./  E0(z)  If  i;0(z)-O 

1.2.  nie  gPHT  for  the  PolBeon  Slatrlbution. 

Suppose  x^,  Xg  ...  Is  a  sequence  of  Independent  ohservatlons 
from  a  Poisson  distribution  vlth  mean  We  want  to  test  the  null 
hypothesis  against  the  alternative  H2:^B^g(>A^)  with  strength 

For  any  positive  Integer  m,  the  Joint  density  function  of  X 

Oi 


under  Is 

m 

(1.2.1) 

•’dA.  ja 

1<b1 

Hence, 

Ex.  -m(A_-A.) 

(1.2.2) 

and 

(1.20) 

z^  «  log  f(Xj^;A2)/f(Xj^;\j^)  »x^  log(A2/^)-(  V^l^ 

Using  Wald's  values  (1.1.4)  for  A  and  B  and  taking  logarithms 
of  (1.1.2)  and  (1.1. j),  the  accepteuice  rules  are  given  as  follows: 
At  the  im-th  stage  of  the  e:qperlment, 

(1)  stop  sampling  and  accept  If 
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m 

(1.2.4)  log  <  log  B* 


11)  atop  sanpllng  and  accept  If 
m 

(1.2.5)  log  (Xg/^)  -  “  ^i)  >  los  A' 


111)  othend.se,  continue  e:g?erlnentatlon. 
(1.2.4)  and  (1.2.5)  are  conveniently  written  as 

log  |o(g/(l-aj^)|  ^  (Xg-Xj^)*! 


(1.2.6)  Z  X  < 
1-1 


log  -  log  X^  log  Xg  -  log  x^ 


n 


s  +  sms  a^^  (say) 


(1.2.7)  2  X  > 
1-1 


log 


I  (l-ag)/aj^ 


(Xo  “  X.,)  m 

+  — = - = - 

log  Xg  -  log  Xj^  log  Xg  -  log 


s  Cg  +  8  m  5  (say)  . 

Thus,  It  Is  seen  that  the  acceptance  lines  are  parallel  straight 
lines  vlth  the  same  slope  s  -  (Xg  -  Xj^)/(log  Xg  -  log  X^^)  vhlch  Is 
Independent  of  the  desired  error  probabilities.  However,  c^  and  Cg 
are  detezxilned  by  the  hypotheses  values  and  the  error  probabilities. 

If  0(^  -  -  a,  from  (1.2.6)  and  (1.2.7)  It  follows  iiaaedlately  that 

Cg  m  -c^m  c  (say). 

1.5  PC  and  ASH  Ptanctlons  for  the  Poisson  3PBT. 

1)  To  obtain  the  OC  function  of  the  Poisson  SPRT  we  have  to  find 
that  value  of  h(x)  (see  (I.I.7)  Tor  which 
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(1.3.1)  2  ^(XgA,)*  e  e“V  x!  »  1- 

aw)  ” 

Sunning  the  series  and  talcing  logarithms,  (I.3.I)  leads  to 


(1.3.2)  -  (Xg-X^)  h  (\)  -  \  -  0. 

Then,  from  (I.I.6),  the  OC  function  L(x)  Is  approximately  given 


by 

(1.3.3) 


However,  It  Is  Impossible  to  explicitly  solve  (1.3. 2)  for  h.  The 
Statistical  Research  Group,  Columbia  IMlverslty  an  In¬ 

direct  method  of  obtaining  a  solution,  but  even  this  method  does  not 
yield  the  OC  function  directly  for  a  given  value  of  \. 

For  practical  purposes,  ve  can  find  the  approximate  values  of 
the  OC  f*inetlon  for  the  following  particular  values  of  X  to  give 
a  rough  picture  of  the  OC  function  of  the  test. 


0  1 


(1.3.5)  X^  1  - 

(1.3.6)  8 

(1.3.7)  X2  0^ 


C^)  m  1/2  If  0^  -  0^  -  O 


(1.3.8)  «  O 

Since  for  X  ••  s  we  have  h(x)  ■  0  as  the  (»Iy  solution  of  (1.3. 
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("by  L'Hospital’s  rule)  (I.3.6)  Is  obtained  by  taking  the  limit  value 
of  (1.1.6). 

17e  shall  introduce  a  grai^cal  solution  for  (l.3>2).  ^  (l.3«S) 


h  log(X.2A3^) 

(1.3.9)  e  ^  ^  -  1  +  (Xg  -  • 

To  solve  (1.3.9)  ior  h,  malce  the  following  transfonoatlons: 


h  log  (X^A,) 

(1.3.10)  yj^  -  e  ^  ,  yg  -  1  +  (kg  -  >^1)  Vk  • 

Since  los{\^/\^)  and  (kg  -  X^)A  are  fixed  numbers  for  any  given 
Xj  y^  is  an  exponential  curve,  and  yg  is  a  straight  line.  Then 
the  solution  of  (1.3. 2)  can  be  found  as  the  intercept  of  the  straight 
line  y^  and  the  exponential  curve  yg  as  shown  below  . 

(k  «*  o) 


The  accuracy  of  the  graphical  solution  can  be  Improved,  if  desired, 
by  a  Newton  iterative  procedure. 

Fraa.  the  h  found  by  this  method  we  can  obtain  the  approxima¬ 
tion  of  L(k)  by  substituting  h  into  (I.I.6). 

11)  In  order  to  obtain  the  ASN  function  free  (I.I.8)  or  (I.I.9) 
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ve  have  to  Imov  the  CX!  function/  L(x).  However/  as  ve  have  seen 
in  the  preceding  parts/  the  OC  function  can  only  he  found  approxl- 
ciately.  For  the  following  particular  values/  the  foUowlnG  approxl- 
natlons  to  the  ASH  function  nay  he  fouz^  from  (1.1.8)/  (I.I.9)  and 

-Cj^/S 

I  (1  -  Oj^)  C^  +  Oj^Cgj.  /(x^-  s) 

■  (1  -  ao)c/(8  -  Xj^)  If  ctj^  ■  0^  ■  o 
-  ■  c®/s  if  0^  ■  0^  ■  a 

I  (1  -  O^)  Cg  +  Q^cJ  /(Xg  -  s) 

-  (1  -  2C»)c/(Xg  -  s)  If  CXj^-Og"® 
(1.3.13)  is  obtained  fron  (1.1.9)  since  E_(*)  -  0.  Iliese  particular 

0 

five  values  of  X  will  serve  sufficiently  veil  for  most  practical 
puxposes. 

2 

Although  it  Is  knovD  that  )  exists/  neither  exact  nor  simple 
approximate  fonoulae  to  obtain  the  variance  of  H  have  yet  been  found/ 
and  very  few  enpirlcal  studies  have  been  undertaken.  Hdwever/  from 
Tables  VII  and  VIII  It  vUl  he  seen  that  the  standard  deviation  of  N 
for  a  truncated  SPRI  tends  to  he  very  large  when  X  Is  between  X^ 
and  Xg  and  a  Is  very  small. 


(1.3. ^^)  to  (1.3.8): 

X 

(1.3.11) 

0 

(1.3.12) 

(1.3.13)  a 

(1.3.14)  Xg 

(1.3.15) 
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1.4.  The  Truncated  SPBT 

Since  the  SPBT  has  no  upper  bound  for  the  sample  eize,  for  prac¬ 
tical  purposer  ve  may  have  to  force  the  test  to  terminate  at  some 
definite  stage,  say  n^,  of  the  e^qperiaent.  Vlald^u7  gave  these 
rules  for  truncation. 

If  the  test  does  not  terminate  by  the  n^-th  stage,  at  the  n^-th 
stage, 

I)  stop  sampling  and  accept  if 

no-1  no 

(1.4.1)  log  B  <  £  >4  ^  1%  A  and  log  B  <  £  s.  <  0 

1-1  ^  1-1  ^ 

II )  stop  sampling  and  accept  Hg  if 

no”!  n^j 

(1.4.2)  log  B  <  £  Z4  <  log  A  and  0  <  £  <  log  A 

1-1  1-1 

In  this  thesis  ve  use  the  following  rules  for  truncating  the 
SPBT.  Let  n^  be  the  maximum  sacgple  sise  of  the  14PRT  (sis  Chapter 
II)  under  the  same  hypothesis  and  strength,  then  the  terminal  decision 
at  the  n^-th  stage  will  be  the  following: 


I)  St<qp  sampling  and  accept  if 

n.  a  4-  b 

£  X.  <  - 2 - ! 

1-1  ^  2 

II)  Stop  sampling  and  accept  IT 


no 

£  X,  > 

1-1  ^  " 


2 
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If  the  SPRT  Is  truncated  at  a  sufficiently  large  n^..  it  Is 
to  he  expected  that  the  error  prohabilitles  vlU  not  he  greatly 
affected.  However «  truncation  will  reduce  the  average  saople  site 
somewhat,  especially  for  Intermediate  \  values  and  it  will  have 
a  similar  effect  on  the  SI3H  function. 


CHAPTER  II 


immiUM  PROBABILITY  RATIO  TEST  AIID  ITS  APPLICATION 
TO  THE  POISSON  DISTRIBUTION 


The  expected  value  of  the  sample  slse  (ASN)  of  a  sequential 
test  depends  on  both  the  paiaoeter  point  0  and  the  particular 
sequential  test.  Ideally^  ve  vould  like  to  find  a  sequential  test 
vhlch  uinlmlses  the  ASIT  functlni  for  all  valoies  of  0>  but  no  such 
"unlfoniily  most  econoDiical"  test  exists.  Hoeffding  given  a 

lower  bound  on  the  ASIT  at  inteniedlate  values  (Chapter  III)  for 
any  sequential  test  meeting  specified  bounds  on  the  error  probablll> 
ties.  In  order  to  come  close  to  achieving  this  lower  bound  under  cer> 
tain  conditions^  Hall  introduced  the  sequential  mlnliaum  proba> 

bllity  ratio  test  (llPRT).  In  this  chapter  ve  give  a  general  dis¬ 
cussion  of  this  test  and  then  consider  the  special  case  of  the 
Poisson  distribution. 

2.1.  The  Miniiaun  Probability  Ratio  Test  (MPHT)  for  Testing  a  Simple 
Hypothesis  against  a  Siqple  Altemative. 

Svqn>086  Xg  . . .  is  a  sequence  of  independent  dbservatlons 
with  a  cccmon  density  function  f(z:0)  and  consider  testing  the  null 
hypothesis  H^:  0  ■>  0^  against  the  altemative  0  ■  Og.  Let  0^ 

be  a  parameter  point  between  0^  and  Og  and  f^  «  f(x:0^)(ia0.1.2). 
Denote  by  3^  «  Xg  ...  x^^)  a  sample  of  size  m  and  f^  the 
joint  density  function  of  X^. 
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We  Introduce  the  weight  functions  and  kg  where  k^  +  kg« 
for  f^  and  fg>  respectively^  to  obtain  a  specified  wel£^ted  average, 
a,  of  the  error  probabilities  (see  section  2.2). 

The  IIFRT  procedure  Is  described  as  follows:  at  the  u-th  stage 
of  the  e^qperlnent, 

I)  stop  saiapllng  and  accept  If 

(2.1.1)  <  a  and  V2»/^l'ln.  5  " 

II)  stop  sampling  and  accept  Hg  If 

(2.1.2)  kiV'on  <  “  VlnA2^2i.  ^  ^ 


111)  otherwise,  continue  sampling. 


Since  the  acceptance  lines  generally  converge  at  some  stage  of 
the  experiment,  say  n^,  the  test  has  the  upper  tound  n^  to  the 
uaxlnum  number  of  trials,  aM  therefore,  either  category  1)  or  11) 
will  occur  for  some  m  <  n.. 


2.2.  Ilajor  Property  of  the  MPRT 

Let  probability  of  making  the  decision 

■  ^2)  when  H^(l  ■  1^2  and  1  J)  Is  true.  Then, 

(2.2.1)  1-.,  oi  .  P,  (a^)  . 

s_ 


1^2  a' 


V  p 
•2  2 


(2.2.2) 


(di)  -  kj 


l6 


vhere  (l  ^  1.2)  is  the  subsets  of  the  sample  space  of  for 
making  decision  by  (2.1.1)  and  (2.1.2)^  respectively. 

Vfe  show  one  of  the  important  properties  of  the  MPBT  by  the 
following  lemma t 
Lenina  1 

For  the  MPRT  with  specified  the  specified  weighted  average 
of  the  true  error  probabilities  is  not  greater  than  the  preassigned 
value  a. 


c 

Proof:  In  S"  it  is  obvious  that 

■  ■'  '  U 


min 

1«1.2 


k.  f . 

1  Im 


'‘l  ^ 


that  is. 


(2.2.5)  >=1  «1  <  «  7  /  “  fo  t''2> 

m 

where  is  the  probability  of  d^^  when  0^  is  the  true 

pareuneter  value.  Similarly,  in 


(2.2.4)  k^  0^  <  P^(d^) 

Renee,  it  follows  immediately  that 


(2.2.5)  ® 

assuming  the  acceptance  lines  converge. 

For  the  case  k^^  ■  kg  ■  1/2,  from  (2.2.5)  to  (2.2.5)  ve  have 
Qtjl^  <  2  a  <  2  a  wid  exj^  +  <  2  a.  If  0^  can 
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"be  80  chosen  that  aj^  <  a  and  Og  <  a. 

In  any  case,  at  least  one  of  the  Inequalities  must  hold  since 
oj^  +  <  2  a. 

Ho  method  of  choosing  0^  to  achieve  specified  Individual  error 
probabilities  Is  known^  except  for  a  limited  result  In  the  case  of 
symmetry  Such  symuetiy  does  not  obtain  In  the  Poisson  case 

considered  here.  We  shall  attempt  to  achieve  equal  error  probabilities 
by  a  method  described  In  section  2.4. 


2.3*  The  MPRI  for  the  Poisson  Distribution 


Let 


be  a  sequence  of  Independent  observations  froa 


a  Poisson  distribution  with  mean  \.  We  vlsh  to  test  the  null  hypo¬ 
thesis  \  against  the  alternative  X  «  Xg  (  >  X^)>  and 

assume  weight  functions  ■  kg  «  l/2  and  a  <  1/2. 

The  MPBT  Is  carried  out  as  follows: 
at  the  m-th  stage  of  the  experiment > 

i)  stop  sempllng  and  accept  If 


u 

It 

1-1 


e  Xj^  I 


m 

jt 

1-1 


~x!  ^ 

V  « 


<  2  a 


l.e. 


(2.5.1) 


a  log  20  m(Xp-  X  ) 

E  — _ _  +  - _ 

!■!  *  log  Xg-  log  X^  log  Xg-log  X^ 


which  Implies  fg^  <  f^^ 

11)  stop  sampling  and  accept  ^  If 


s  c^+  r^m  (say) 


-  log  20! 


(2.3.2)  2  X,  > - 

1.1  ^  “  log  -  log 


•f 


log  log 
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S  C2  +  TgH  (say) 


vhich  Implies  >  f 

111)  otherwise,  continue  eanpllng. 

(2.3.1)  and  (2.3.2)  will  define  the  acceptance  lines  of  the  test. 
These  two  lines  meet  at  m  .  n*  where 


- 


^2 "  ^1 


and  hence,  this  Is  the  upper  hound  for  the  sample  size.  !nie  actual 

maximum  sangple  nvimber  n^  will  be  an  Integer  sllf^tly  smaller  than 

n'. 

o 

An  Illustrative  diagram  for  carrying  out  this  test  with  a  coo^ 
parlson  of  the  truncated  SPRT  and  the  FSST  appears  in  Figure  1 
In  the 


2.4.  A  Choice  for  0 

o 

In  the  I4PBT  It  Is  possible  to  choose  0^  at  any  parameter 
point  other  than  0^  and  Og,  depending  on  our  emphasis  on  the  error 
probabilities  and  and  the  vei£^t  functions  and  kg. 

In  this  thesis,  ve  will  choose  In  an  attempt  to  obtain  approxi¬ 
mately  equal  error  probabilities  at  0^,  l.e.  Pq(^)  ■ 
this  point  of  view,  ve  take  0^  as  the  value  which  yields  equal  "di¬ 
vergence"  /J7  between  0^^  and  ©^,  and  between  0^  and  Og  . 

Suppose  the  null  hypothesis  is  }I^:  f(x)  .  f^(x)  and  the  alter¬ 
native  I^:  f(x)  s  f2(x)/  ^or  a  given  x,  the  mean  information 
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for  discrir-ilnatlon  per  obaervation  in  favor  of  acainst  Hg  ia 
aef  ined  £'(]  as 

p  f,(x) 

I  (l;2)  ■  /  f,(x)  log  -  d  X 

^  ^  fg(x) 


^ere  K  ia  the  entire  aaiqple  qpaee.  I  (2:1)  ia  similarly  defined. 
Under  the  assumption  that  the  probability  measures  for  f^(x)  and 
f2(x)  are  "absolutely  continuous"*  the  integrals  always  exist. 

The  divergence  D(1«2)  between  and  is  defined  /  jy  as 

f  Mx) 

(2.4.1)  D(l,2)  .  I(l:2)  ♦  I(2tl)  -  /  JIM  -  f«(x)  7log  d  x. 

^  ^  ^  fg^x) 

D(1*2)  is  a  measure  of  the  difficulty  of  discriminating  between 
hypotheses  and  Kg*  and  thus  will  be  a  reasonable  choice  to  deter¬ 
mine  0  for  the  MPRT.  Further  detailed  discussion  concerning  in- 
0 

formation  and  divergence  will  be  found  in  KuUbadc  / 


Exaiqple:  Poisson  case 


Suppose 

X 

-  X^ 

i.e. 

f^Cx) 

•  xj^  e  Vxl 

and 

Hg:  X 

l.e. 

fgC*) 

.  •  \x'. 

Then*  by  definition 

m 

A. 

.’"i 

xj  .'Nx! 

I(l:2) 

-  z 

XeO 

xl 

log 

M 

•K 

Z  -= 
XBO  x'. 


e 


x(log  -  log  Xg)  -  (Xj^  -  Xg) 
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(2.4.2)  a  (log  -  log  Xg)  Xj^  -  ()^  -  Xj) 
SlinilBrlyi 

(2.4.3)  1(2:1)  =  (log  Xg  -  log  X^)  Xg  -  (Xg  -  \^) 
Hence,  from  (2.4.1)  to  (2.4.3), 


(2.4.4)  D(l,2)  -  (Xg  -  Xj^)(log  Xg  -  log  \^) 

For  some  X^,  where  ^  ^  difficulty  of  discrlni- 

natlng  between  hypotheses  X^  and  X^  and  between  X^  and  Xg  is 
equal.  That  is,  X^  is  obtained  from 
D  (1,0)  -  D(0,2)i  that  is 

X-Q  -  log  Xj^)  »  (Xg  -  XQ)(log  Xg  -  log  Xjj) 


(2.4.5)  g(XQ)B  (log  Xg-log  >^)Xq  +  (Xg-Xj^)log  Xjj  +  (X^log  X^-Xglog  Xg) 


m  0  . 


Newton's  approximation  method  was  used  to  solve  this  equation 
numerically,  starting  with  the  initial  value  ^  *  (Xj^  +  Xg)/2. 
That  is. 


^0.1+1 

(2.4.6) 


■  X 


0.1 


o 

X.o(lo6  V  N^-X-glog  Xg 

log  Xg  -  log  Xj^  +  (Xg  -  \)/\ 


(1  -  0.1.  ...) 


with  the  criterion  for  stopping  the  Iteration 


I  ■  \.i  I  5  “ 


,-5 
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and  i  <  200. 

Sy  thlB  method,  divergence  values  In  the  Appendix  vere  conq^ter 

calculated.  Thla  X.  value  is  there  denoted  hy  D. 
o 

A  few  exBEiples  were  considered  using  ®  instead  of  «  D; 

a  cottqparlson  of  these  choices  is  discussed  in  Chapter  V.  It  was  found 


that  D  <  8. 


CHAPTER  III 


FIXED  SAtIPLB  SIZE  TEST  AND 

HOEFPDING'S  LCVIER  BOUND  FOR  THE  EXPECTED  SAIIPLE  SIZE 


5.1.  Moat  Powerful  Test  for  t!^ Pel8ttoti:iPl0tglb6tlon. 

To  laake  a  saatple  size  conparison  between  the  MPRT  and  a  fixed 
sanqple  size  test  (FSST)^  ve  need  to  find  the  mlnlTmim  aangple  size  for 
the  Eiost  powerful  FSST  with  specified  hounds  on  the  error  probabili¬ 
ties.  The  most  powerful  FSST  Is  obtained  by  the  Heymn-Pearson 


leimia. 


For  the  Poisson  case  with  sangple  size  m>  the  test  Is  carried 
out  as  follows: 


1)  accept  Hg  If 


f^  ^2  ®  '*1  * 

(3.1.1)  - - -  >  k(>0) 

f,  .xi  ^1/  , 

lu  n  X,  e  /x  I 

1-1  ^  ^ 

11)  otherwise j  accept  H^^  . 


Ibking  losarlthos  on  both  sides  of  (3.1.1)^  It  follows  that 


m 

Xj  (log  Xg-  log  Xj^)  -  n(Xg  -  Xj^)  >  k' 


(3.1.2)  Z  X,  >  k" 

1-1  ^  “ 

Q 

We  know  that  the  randoa  variable  y  «  Z  x,  has  the  Poisson  dis- 

1-1  ^ 

tributlon  with  mean  iaX>  l.e. 
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(5-1.5)  s(y)  ■  (“  e"“^/y  ! 

If  the  Type  I  error  and  the  size  of  san^jle  m  are  given, 
k"  Is  obtained  as  the  suallest  Integer  for  vhlch 

(5-1.4)  1  -  (k-  -  1)  < 

where  F^(^)  is  the  Poisson  distribution  function  at  x.  If  the 
Type  I  error  and  the  Type  II  error  are  preassigned,  the 
sas^le  size  can  be  obtained  "by  finding  k"  and  the  minimal  m  vhlch 
will  satisfy, 

(5.1.5)  1  -  (k"  -  1)  <  otj^ 

and 

(M.6)  (k"  -  1)  < 

5-2.  Normal  Approximation  for  the  Poisson  FS3T 

If  \  Is  large,  the  Poisson  distribution  can  be  approxliaated  by 
the  iTonoal  distribution.  That  Is, 

a  «  (^  ~  ^  V?  ) 
t 

-  f  (1/  /Si  ) 

Then,  by  (5.1.5)  and  (5.1.6) 

(5.2.2)  9  ( - i - ) 

/5^ 


-x®/2 

e  '  d  X 


>  1  -  Otj^ 


(5-2.1)  F^ix) 

where 

♦  (t) 
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k”  -  1  -  mX  +  1/2 

(5.2.5)  9  ( - — -S - )  <  OL 

are  the  required  Inequalltlee. 

If  the  Poisson  mean  Is  not  larger  the  normal  approximation  may 

1/2 

not  be  adequate  so  that  the  following  rule  ma  used.  If  (m  X^)  ' 

<  4.0j  several  integer  values  were  chosen  in  the  neighborhood  of 
the  m  \dilch  was  obtained  from  the  Normal  approximation  method. 

These  were  substituted  into  (3.1.?)  and  (3>1>^)>  vlth  suitably  chosen 
k"  s  from  the  Poisson  tables^  until  (3.I.?)  and  (3.1<6)  were  satis¬ 
fied.  The  minimum  such  m  is  the  sample  size  of  the  most  powerful 

test.  In  the  tables  In  Appendix  those  values  obtained  by  the  Hozual 

#  # 

approximation  appear  with  such  as  n^  and  FS9T  . 

3.3*  Hoeffdlng's  Lower  Bound  (Hlfi)  /~27 

Let  x^,  Xg  . . .  be  a  sequence  of  Independent  observatloos  with 
a  common  density  function  f(x:0),  and  suppose  we  want  to  test  the 
null  hypothesis  0  •  against  the  alternative  Hgt  0  ■  Og. 

Also,  let  f(x:0)  »  fj^(x;0)  under  Hj^(l  ■  1.2),  and  the  decision 
Is  to  accept  Hj^. 

As  one  of  the  optimum  properties  of  the  SPBT,  It  has  been  proved 
^1£7  'the  SPRT  minimizes  the  expected  sample  size  at  the  points 

0^  and  Og  subject  to  specified  bounds  on  the  error  probabilities. 

In  general,  Its  expected  sample  size  Is  largest  when  0  Is  between 
0^  and  Og  but  not  always  (see  Figure  2).  However,  there  exist 
tests  \dilch  have  a  smaller  expected  saiple  size  at  intermediate  0 
values  than  the  SPBT.  Kiefer  and  l/eiss  have  discussed  isipor- 


tant  qualitative  properties  of  such  tests.  These  tests  loay  be  Judged 
by  con^rlng^  at  any  parameter  point  9,  the  expected  s:\nple  size  of 
the  test  with  the  smallest  expected  sample  size  attainable  by  any 
test  having  the  same  error  probabilities  at  0^  and  O^. 

The  Hoeffding  Lower  Bound  on  the  ejg^cted  sample  size  at  any 
0^  is  given  by 

(3.3.1)  Eg  (N)  >  Slog 

assuming  that  the  following  integrals  exist 


+  -r/kY/i 


(3.3.2)  s  -  (Si^Cg)  >  ^1  “  /  ^o  ^V^i^  ^ 

and 

(3.3.3)  -  /{  loe  (fg/^i)  -  ^1  +  ^2  'o  • 

Also  assume  that 


f  (x)  a  0  Implies  min  f,  ■  0 
°  i»1^2  ^ 


and  that 

Hop 
B-  (B  Y.r  “  T  E  (H) 

J  r  1 

Is  satisfied  where  Yj  ■  log  |  f2(x)/fj^(x)  |  "  ^1  ^2 

Hoeffding  proved that  f or  1 11  (3.3.1)  gives  a  lower 

bound  among  all  strength  (o^.o^)  tests. 

For  the  Hoxmal  density  function  with  variance  one  and  mean  0^ 
where  0^  *0^  0^  >  -&  and  0^  »  5  >  0,  be  coopered  the  numerical 
values  obtained  by  his  lower  bound  with  those  of  the  fixed  8eBg>le 
size  test,  one  of  Anderson's  tests  /”l7  'the  SPRT  of  the  same 

strength  (o^^,  o^) .  For  0;^  ■  0^  a  a  <  1/2,  the  results  Indicate  that 
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equality  in  la  nearly  obtainable  with  a  FSST  if  a  la  very 
email  and  vlth  the  SPRT  If  a  la  auffielently  large.  For  other  a 
valuea^  thoae  cooiaonly  uaed  In  practice,  the  IIPRT  (vhlch  colncidea 
vlth  one  of  Anderaon'a  teata)  nearly  attaina  thia  bound. 


Hoeffdlng'a  Lower  Bound  for  the  Polaaoo  Distribution 


Suppoae  f(x:\)  >  \  e  /x  !  and  ve  wish  to  teat  \ 
agalnat  X  ■  Xg  (>  X^)  with  strength  of  the  test  ■  a 

<  1/2.  Then,  by  definition  (5*3*2) 


XbO 


^0  ®  *  log(XQ/Xj^)-(Xjj-Xj^)  E 

XbO 


(J.^^.l)  -  Xq  (log  Xg  -  log  X^)  .  (x^  -  x^). 

Therefore, 

(5.4.2)  Cl  -  (log  X^j  -  log  \^)  -  (Xq  -  Xj^) 


(3.4.3)  Cg  -  X^  (log  X^  -  log  Xg)  -  (Xq  -  Xg) 
Also  by  definition  (3. 3 >3) 


-  (log  log  \)  +  (Xp  -  Xj^) 


+  (log  X^  -  log  Xg)  +  (Xg 
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CO 

a>  S 

x»o 


X 


1 


x(log  X.g-  log  \j^)+  \Q(log 


-  log  \g)  j- 


00 

Z 

XbO 


log  (x 


■  Xq  (log  Xg  -  log 


For  aoy  X^  ve  can  obtain  Boeffdlng's  lower  bound  by  substituting 
(3.^^.2)  or  (3*4.3)  and  (3*4.4)  into  (3*3.1)* 

Froo  (3.4.2)  and  (3*4.3)  it  is  easily  seen  that  the  relation 
■  ^2  a  ^  is  obtained  at 


X^j  -  (X^  -  X3^)/(log  Xg  -  log  \^) 


and  this  value  of  X^  Is  equal  to  b,  the  comon  slope  of  the  accept¬ 
ance  lines  of  the  SPRT  for  the  Poisson  distribution  under  the  same 


hypothesis.  For  X  ■  s,  Hoeffding's  lower  bound  Is  given  by 

fr  ^  ^ 

-  (s  log  ^  -  «-ii)iog  a*]- 


■11/2  1/2 
■T  log 


(3*4.5)  E  (R)  > 


S  log  Y  -  8 


The  E  (H)  of  the  SPRT  Is  given  approxinately  by  (1.3*13)2  therefore, 

0 

we  can  coopare  Hoeffding's  lover  bouzid  with  the  approxiioatlon  of  the 


SPRT  ASR  at  X^  ■  s. 

o 

The  HIB's  In  the  tables  In  the  Appendix  were  cooqputer  calculated. 


CHAPTER  IV 


THE  PROGRAMS  FOR  THE  I4PRT  AND  THE  SPRT^ 

0 

FOR  THE  POISSOI^  OiaTRIBUTION 

In  thlB  Chapter  ve  dlsciiss  the  proszams  by  vhlch  the  exact 
OC,  ASH  and  SM  functions  of  the  MPRT  and  SPRT^  vere  obtained  for 
the  Poisson  distribution.  The  prograus  (for  the  IIPRT  and  the  SPRT^) 
were  written  In  the  "IT”  langviage  and  stored  as  "K.  Fukuahliaa, 

MPRT  -  a"  and  "K.  Kukvishlna,  SPRT  -  A">  respectively,  In  the  library 
of  the  Reseazvh  Coc^tatlon  Center,  the  consolidated  University  of 
North  Carolina,  Chai>el  Hill,  North  Carolina,  for  future  use.  Ihe 
detailed  coqpller  prograu  for  the  MPRT  Is  shown  In  A.?. 

4.1.  Brief  Explanation  of  the  Prograia  for  the  MPRT 
Let 

Bij  (integer):  the  acceptance  boundary  for  the  1-th  trial. 

(integer);  the  rejection  boundary  for  the  1-th  trial. 

p(q^):  the  probablll^  of  acceptance  at  the  1-th  trial. 

p(^):  the  probability  of  rejection  at  the  1-th  trial. 

p(m.  the  probability  of  J  defects  (j  »  E  x)  at  the  1-th  trial. 

1.  j 

n^:  the  point  at  which  the  acceptance  and  rejection  lines  Intersect, 
n^:  the  naxliiiuD  possible  number  of  trials;  the  least  1  such  that 

“l  “  “i  5  1 
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n^:  the  number  of  Poisson  probabilities  to  be  calculated. 

Hi  «  Cg  -  Cl  +  5 

Sgi  the  ox^nates  (L  x)  of  the  acceptance  and  rejection  lines 

at  n'  . 
o 

X  1,  \  2>  ...XU  (input):  The  parameter  points  for  iihleh  the  0C« 

ASS  and  SDN  functions  are  to  be  calculated. 


By  (1.2.6)  and  (l.2.7)>  after  ve  detexmlne  the  boundaries  for 
acceptance  and  rejection  and  n^>  ve  perform  the  foUoving  calcula¬ 
tions  for  (i  <  n^) 

—  o 

ialn,(j,i^  j^-1) 

vfaere  F^Cx)  is  the  Poisson  probability  of  x  vlth  mean  X. 


P(m^) 


m  -k 

{"'(“n.i)-  } 


k«n.  ,+l 

—V*  1 


PAtfg 

for 


A  CCMPAEaSQH  OF  SEqjJENTIAL  TESTS  FOR  THE  POISSCN  PARAMEZEER 

tjy 

Kozo  Fiikusbina 


Page 

7 

15 


line  tnna- 
ber  from 

the  top  Misprint 

7  AIqCz)  if  EqCz)  =  0 


Correctlcn 


/Bq(z^)  If  E(z)  *  0 


17 


17 


$  i  1  ^  2  s 

i5i  \ 


7^ 


/x'l 


fo-/f  * 
2a'  cm 


ill  >^2  ®  /*i' 


m 

l5l^ 


-K. 


Vx^S 


17 

20 

vAiich  teiplles  f^^^  <  f^ 

vhich  impUeB  <  f 

19 

16 

i.e.  f2(x)-\^e  ^/xl 

r  “^0 

i..eo  fg(x)  «  \g  e  /xl 

20 

16 

”  '‘o.i  •  IT 

0 

^o.i+1  ‘  ^Ooi  “  "ST 

o 

20 

17 

i2vk^6)  ''■■  + 

(2.4c6)  r:  - 

33 

20 

The  tlae  for  conpitatlcn 

Obe  time  for  coD^Uatlon  I 

36 

16 

the  MPRT  Is  imifOTnOy 

Che  SlIT  of  the  MPFT  is  unifoxmly 

33 

5 

(”os  -  ”) 

36 

2 

(n«,  “  5T1) 

(%S  ”  ^  J 

67 

17 

.oooTn  0005',  .00025 

.0000071  .000034  .000023 

l\ 
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Tlie  overall  flow  dlegraa  Is  shown  on  the  following  page. 
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i 


4.  The  steps  A/  B  •••  M  appear  In  the  compiler  prograia  in 


A.  7  correspondingly. 
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4.2.  Variable  Asalgntaents  for  «  A" 

Y  1  (input)  m  Y  2  (input)  -  Y  5  (input)  »  Xg 

Y  4  (input)  -  a,  Y  28,  Y  29  ...  Y  38  (input)  ■  X  1,  \2,  ...  X  11 
N  1  -  i,  N  2  -  3,  N  900  -  n^,  N  902  »  n^ 

N  10  —  N  500  I  ,  N  500  -  H  904  t 

N  905  —  K  1050  :  Alplianuaeric 

Y5-Cj^,  Y6«rj^,  Y7-Cg,  YS-r^ 

Y  100  —  Y  499:  #  y  500  —  Y  999  : 

Y  1000  —  Y  1999^  Pj^(j)  *  Y  2000  -  Y  2400:  P(m^) 

*  i 

z  0  —  z  999  *  E  p,(x)  >  z  1000  —  z  1999*  z  P.(x) 
xmi  ^  x»o  ^ 

Z  2000  —  Z  2400  t  P(mj^)  vhere  J  «  0,  1,  2,  . . .  n^^ 

4.3.  The  Proeram  for  "SPOT  -  A" 

The  storage  spaces  for  the  variable  assignments  are  very  siml> 
lar  to  the  program  for  the  I<!PBT,  except  for  the  following  changes. 

For  input, 

i)  Yl-Xj^,  Y2-X2,  Y3-o^,y4-q^ 

11)  the  truncation  point  (H  90O  ■  n^)  must  be  given 
ill)  the  noober  of  paraxieter  points  for  which  the  OC,  ASH  and 
SM  functious  are  to  be  calculated  is  10(Y29,  • . .  Y38)  Instead  of 
11  in  "MPRT  -  A"  . 

4.5.  Outputs  and  Capacities  of  the  Programs,  *YlPRT  -  A"  and  ”3PRI»A". 
unconditional  outputs  for  these  programs  are  as  follows: 


33 

X-Q#  Xgi  0!  &,nd  Xl#  \2  •  •  •  \11 

’'^1'  °2*  ^2*  “o 

c^,  s  and  Cg  (for  the  SPRT  -  A) 
a^  and  ag  (where  a^^  ■  ag) 


For  X  » 

XI,  X2,  ...  X  11 

n 

n 

0 

0 

E 

1-1 

\(n  )  ,  E 

^  ^  1»1 

ASN>  the  variance  of  N>  and  SEN 

For  conditional  outputs  of  the  procrams,  the  following  outputs  are 
added: 

n^  ,  and  (i  »  0,  1,  2,  ...  n^) 

J  « 

>  2  P  (x)  and  E  P  (x)  where  -O,  1,  2,  ...  n. 

xao  x-J 

P(ia.  .)  for  the  1-th  trial  (l  «  1.2  ...  n  ) 

1.  J  0 

P(n  )  and  P(o  )  (for  SPRT  -  A) 

““o  “o 

By  "IIPRT-A"  and  "SPRT  -  Pi*,  we  can  take  any  hypothesis  values 
end  Xg  and  the  error  prohabllltles  a  <  -g  as  long  as  n^  and 
E  X  do  not  exceed  400  and  ^00,  respectively.  Also*  the  core 
storage  used  In  the  computer  Is  7^  for  "MPBT  -  A"  and  7752  for 
"SPRT  -  a".  However^  one  can  change  the  program  according  to  the 
reci:ulrement8  of  each  particular  problem  and  a  computer  capacl*^. 

The  tine  for  cocputatlcn  was  about  three  minutes^  and  the  calcu¬ 
lation  of  the  0C>  ASN  and  SDN  functions  at  one  parameter  point  for 


each  of  eiglit  sets  of  hypotheses  vlth  fovir  error  levels;  tliat  ls> 
one  \  value  for  each  of  thirty- tvo  test  situations,  todc  about 
fourteen  minutes.  Similar  time  Is  required  for  "SPRT-A". 


CHAPTER  V 


SUMIIARY 

In  this  chapter  ve  discuss  the  characteristics  of  the  MPBT 
frou  the  data  obtained*  with  cooparison  of  other  tests*  the  SPBT^* 
the  SPRT  and  the  FSST. 

By  the  use  of  dlvercence  to  obtain  X^*  and  =  kg  ■  l/2,  the 
test  nearly  achieves  ®  small, 

and  and  Xg  are  not  very  distinct.  As  already  mentioned  in 
section  3*2,  the  sum  of  the  exact  error  probabilities,  * 

is  smaller  than  the  preasslcned  level  2a  .  Moreover*  even  thoueh 
the  test  for  the  Poisson  distribution  is  not  symmetric*  we  observe 
tliat  5  0^  5  ot  is  not  very  small*  o'  <a4 

in  General.  As  mentioned  in  the  introduction*  the  use  of  X^  ■  D 
Gives  sliGhtly  better  results  to  achieve  PQ(dj^)  » 
larce  a  values  and  X^  ■  S  g1v®s  sliGhtly  better  results  for 
small  a. 

Cccrparinc  the  OC  functions  of  the  MPRT  with  the  SPHT^  and 
FS9I  the  foUowlnc  points  may  be  found: 

i)  the  SPRT^  has  Generally  bicker  discrimination  than  the  MPRT 
between  X^^  and  Xg  for  small  a*  the  OC  function  of  the  MPRP 
tends  to  be  close  to  the  SPBT^.  However*  the  difference  between 
the  OC  functions  is  not  sufficiently  larce  to  be  of  particular  im¬ 
portance*  and 

li)  the  l^tPHT  has  Generally  hiGher  discrimination  fOr  small  a  than 
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the  PSST  except  near  and  Xg.  But  if  a  is  larce^  the  PSST 
tends  to  have  uniformly  hl^er  discrimination. 

From  the  tables  aud  shorm  in  the  Appendix  ve  see  the 

follovlnG  characteristics  of  the  ASN  function  of  the  14FBT. 

I)  Tlie  smaller  the  Chvalue>  the  closer  to  s  Is  the  iiaxiuum  value 
of  the  ASN  of  the  MPRT. 

II)  For  small  a,  the  ASH  of  the  MPRT  Is  smaller  than  both  the  SPRT^ 
and  SPRT  for  some  values  of  X  between  X^  and  Xg>  but  the  maxinun 

of  the  ASH  for  the  MPRT  and  SPHT^  are  not  very  different, 
ill)  The  ASN  of  the  MPRT  is  unifonaly  smaller  than  for  the  PSST  ex¬ 
cept  for  extremely  small  a  values. 

Iv)  The  HLB  at  X  «  D  (and  X  =  •)  is  nearly  attained  for  large 
values  of  a  by  the  MPRT.  For  smaller  values  of  a>  even  though 
the  ASN  of  the  MPRT  Is  not  close  to  the  HLB«  It  Is  closer  than  the 
ASI?  of  the  SPRT  or  the  PSST. 

We  observe  that  for  any  a  values  the  MPRT  Is  unlfomly  smaller 
than  for  the  SPRT^  which  is  presumably  smaller  than  for  the  SPRT. 

Therefore,  the  MPRT  ui.ler  these  ccndltloos  appears  to  be  ad¬ 
vantageous  as  compered  with  these  other  tests  when  the  paretneter 
point  lies  near  the  average  of  X^  and  Xg,  and  particularly  idien  a 
Is  small.  Ever  If  the  ASN  of  the  liPRT  is  8ll(^tly  larger  than  of  the 
SPRT,  the  use  of  the  MPRT  may  be  reccosended  because  of  its  smaller 
SHI.  Further  studies  of  the  MPRT  for  various  weight  functions 


and  kg  should  yield  useful  results. 
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APPEIIDIX  I 

A.l.  Notation 

The  foUowlnc  notation  is  used  In  the  tables  and  flcures 
in  the  Appendix: 

OC  :  the  pperatlnQ  characteristic  function 

ASn  :  the  averase  saicple  tmciber  function 

SDtT  :  the  standard  deviation  of  the  saisple  size  N 

FSST  ;  the  nost  powerful  fixed  saixple  size  test 

HI£  :  Hoeffdinc's  lower  bound  for  the  ASN 

MP6T  :  the  mlnimun  probability  ratio  test  ■  D) 

the  nunbers  in  parentheses  in  the  tables  V  and  VI  were  obtained 
by  the  IIPOT  (X^  -  s) 

SPBT  :  the  sequential  probabilll^  ratio  test 

SPRI  :  the  sequential  probability  ratio  test 

the  numbers  which  have  *  in  the  SPRT.  in  the  tables  VII  and 

o 

VIII  were  obtained  by  hhld's  approxlriatlon  for  the  SPBT 
SPBT^:  the  sequential  probablll'ty  ratio  test  truncated  at  n^ 
a  :  the  specified  bound  on  each  error  probability 

D  :  X^  value  for  which  the  dlversence  between  X^  and  X^  equals 

the  divergence  between  X^  and 
s  :  the  slope  of  the  SPBT  acceptance  lines 

n^  :  the  sample  size  of  the  FSST  obtained  from  the  Poisson  distri¬ 
bution  (with  linear  interpolation  In  the  tables) 
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Up  :  the  sastple  size  of  the  FSST  obtained  froD  Norrial  approximation 

n^  :  the  maxlriun  sanqple  size  of  the  MPBT  (\^  >  D) 

:  the  maximun  sample  size  of  the  t^RT  (\^  ■  :$) 

tip  :  the  HLB  at  X  >  D 


the  wra  at  X  >  s 
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APPENDIX  III 

A. 5*  TABLES  I  -  IV;  Characteristics  of  the  Poisson  MPRT  and 

SPOT  Approximation. 

A. 3.1.  TABIHS  I  -  1,  2,  3,  k. 

X  a  .1  against  HgS  X  ;  .3 

s  «  .18205  D  =  .18652 


TABLE  I  -  1 

a  =  .1,  n^  «  52,  iip  =  31 


X 

OC 

ASN 

HIB 

SDN  of 

MPOT 

MPOT 

SPRP» 

SPOT* 

0 

1. 

1. 

15. 

11. 

0. 

.1 

.9203 

.90 

22.83 

19.75 

7.88 

.17 

.585U 

24.65 

10.37 

s 

.5175 

.50 

24.28 

21.97 

20.40 

10.62 

D 

.U929 

24.11 

19.18 

10.71 

.25 

.2112 

20.51 

10.97 

.3 

.0932 

.10 

16.90 

13.73 

10.13 

TABIE  1-2 
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TABIE  1-4 

a  =  .001,  =.  223,  n#  =  180 


\ 

OC 

ASN 

HIB 

SDN  of 

MPRT 

B 

SPHT» 

SPRP* 

0 

1. 

1. 

55. 

35. 

0. 

.1 

.9991 

.999 

95.04 

76.47 

22.04 

.17 

.6570 

142. 18 

38.03 

8 

.5028 

.50 

l4i.i 

217.1 

130.8 

40.22 

D 

.4459 

139.7 

121.1 

41.13 

.25 

.0256 

92.25 

39.04 

.3 

.00088 

.001 

65.01 

53.19 

27.11 

A.3.2. 


TABLES  II  -  1.  2,  3.  4. 

\  B  .3  against  \  «  .8 
8  ■  .63829  D  B  .64122 


TABIE  II  -  1 


a  -  .1,  -  87  ■  ^7 


\ 

oc 

ASS 

HI£ 

SDN 

MPRT 

MPBT 

SPRT* 

IIPRT 

SPOT* 

0 

1. 

1. 

11. 

8. 

0 

.5 

.9127 

.90 

31.14 

27.38 

13.22 

.6 

.6448 

36.09 

15.40 

a 

.5064 

.50 

36.07 

34.24 

31.56 

15.80 

D 

.4958 

36.03 

30.85 

15.83 

.7 

.2991 

33.84 

15.99 

.8 

.0923 

.10 

27.29 

23.42 

14.53 

TABI£  II  -  2 

a  a  .05#  *  126  Dp  ■  78 


\ 

OC 

ASN 

_ 

HIS 

SEN 

MPRT 

A 

MPRT 

SPRT* 

IIPBT 

SPRT* 

0 

1. 

1. 

15. 

10. 

0 

.5 

.9559 

.95 

46.87 

40.88 

18. 4l 

.6 

.6773 

58.53 

22.51 

a 

.5042 

.50 

58.93 

61.49 

55.15 

23.23 

D 

.4908 

58.85 

51.82 

23.28 

.7 

.2491 

54.14 

23.64 

.8 

.0458 

.05 

40.66 

54.96 

20.32 

TABLE  II  -  5 

a.  »  .01,  «  217  Dy  »  155 


oc 

ASN 

r . .  ■" ' 

RU3 

SDN 

MPRT 

tIPRT 

SPRT* 

■QB' 

1. 

1. 

25. 

16. 

0 

.9911 

.99 

82.02 

69.29 

27.89 

.7368 

116.5 

38.15 

.5018 

.50 

119.3 

149.75 

110.4 

39.50 

.4831 

119.1 

107.3 

39.62 

.1693 

104.8 

40.82 

.0090 

.01 

70.18 

59.26 

30.58 

TABLE  II  >  4 

a  -  .001,  -  550  nj  -  274 


OC 

ASN 

Hlfi 

SDN 

MPRT 

HPRT 

SPRP* 

ij>Rr 

8PRT* 

1. 

1. 

4o. 

24. 

0 

.9991 

.999 

150.4 

106.1 

36.81 

.7973 

205.5 

58.43 

.5007 

.50 

215.5 

338.3 

201.9 

60.06 

.4760 

1 

! 

214.9 

195.8 

60.36 

.1025 

1 

179.3 

63.20 

.00069 

.001 

110.5 

90.69 

40.10 

1^5 


A, 3.3.  TABLES  III  -  1.  2.  3.  4. 

\  »  .5  against  Hgt  \  »  1. 

I  «  .72135  D  -  .72850 


TABLE  III  -  1 

0!  ■  .1#  n^  »«  34  Dp  s  22 


\ 

OC 

AS 

HLB 

1  SDN 

MPRT 

MPRT 

SPRT* 

MPRT 

0 

1. 

1. 

6.0 

5. 

0 

.5 

.9185 

.90 

13.81 

11.60 

5.40 

.65 

.6683 

15.70 

6.45 

s 

.5083 

.50 

15.62 

15.95 

12.87 

6.74 

D 

.4929 

15.58 

12.41 

6.76 

.85 

.2553 

14.11 

6.86 

1. 

.0877 

.10 

11.47 

9.21 

6.22 

TABIE  in  -  2 

a  »  ,05,  Dq  ■  52  nj.  “  52 


\ 

OC 

ASN 

HIB 

SDN 

MPRT 

IIPRT 

SPRT* 

MPRT 

SPRT* 

0 

1. 

1. 

9. 

6. 

0 

.5 

.9594 

.95 

20.33 

17.32 

7.49 

.65 

.7054 

24.89 

9.29 

8 

.5084 

.50 

25.02 

25.02 

21.67 

9.76 

D 

.4884 

24.95 

20.84 

9.80 

.85 

.2026 

21.95 

10.00 

1. 

.0442 

.05 

16.77 

15.76 

8.66 
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TABLE  III  -  5 

a  ■  .01,  n  a  88  >  64 

o  P 


\ 

oc 

ASN 

HI£ 

SDN 

IIPOT 

MPRT 

SPOT* 

MPOT 

SPBT* 

0 

1. 

1. 

15. 

10. 

0 

.5 

•  9915 

.99 

35.1^2 

29.35 

11.19 

.65 

.7655 

48.69 

15.51 

s 

.5014 

.50 

‘^9.85 

60.93 

45.01 

16.48 

D 

.4736 

49.66 

43.11 

16.59 

.85 

.n85 

40.82 

17.04 

1. 

.0084 

.01 

1 

1 

28.22 

23.32 

12.81 

TABLE  in  •  4 

a  -  .001,  »  142  nj.n, 


\ 

OC 

ASN 

HI£ 

SEN 

MPOT 

MPOT 

SPOT* 

MPOT 

SPOT* 

0 

1. 

1. 

23. 

14. 

0 

.5 

.9991 

.999 

56.07 

44.95 

14.71 

.65 

.8289 

84.77 

23.53 

8 

.5006 

.50 

89.25 

137.64 

82.34 

24.92 

D 

.4640 

88.87 

78.61 

25.17 

.85 

.0594 

67.96 

25.76 

1. 

.00082 

.001 

44.07 

35.69 

16.71 
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A.3.4. 


TABIES  IV  -  1.  2.  3.  4 


H^!  \  «  5*  against  Hgt  \  ■  8. 

a  «  6.3829  U  «  6.4122 


I 


TABLE  IV  -  5 


a  B  .01  •  n  B  22f  n_  b  i6 
o  P 


00 

Asn 

MPRT 

SPBT* 

UPRT 

SPOT* 

1. 

1. 

5. 

2. 

.9929 

.99 

8.92 

6.93 

.7447 

12.82 

.5070 

.50 

13.28 

14.98 

.4806 

13.11 

.1605 

11.56 

.0070 

.01 

7.76 

5.93 

11.04 

10.75 


0 

2.94 

4.03 

4.09 

4.19 

4.31= 

3.18 


TABLE  IV  -  4 

a  B  .001,  B  33, 


00 

ASS 

IIPOT 

SPOT* 

MPOT 

SPOT* 

1, 

1. 

4. 

3. 

.9995 

.999 

13.74 

10.60 

.8010 

21.70 

.4989 

•50 

22.80 

33.85 

.4735 

22.75 

.0969 

19.00 

.00069 

.001 

11.81 

9.07 

20.19 

19.58 
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APPETJDIX  IV 

A. 4.  TABLES  V  -  VI;  CharacterlBtlCB  of  the  Pole son  1#RT 

>  D  and  ■)  and  SPBT  Approxlmtlon 

A. 4.1.  TABIBS  V  -  1.  2.  3.  4 

H^3  \  >  .1  agalnat  B^:  \  b  .3 

B  -  .24833  B  -  .26129 


TABI£  V  -  1 

a  «  .1,  n^  «  18,  (n^  .  18),  n^,  -  11 


\ 

oc 

A&tf 

HUB 

UPRT 

SPOT* 

MPOT 

SPOT* 

WPOT 

0 

1. 

(1.) 

1. 

7. 

(7.) 

6. 

0. 

(0.) 

.1 

.9378 

(.9392) 

.90 

9.08 

(8.92) 

7.44 

2.80 

(2.56) 

.2 

.6856 

(.6933) 

9.75 

(9.51) 

3.66 

(3.43) 

8 

.5391 

(.5496) 

.50 

9.50 

(9.28) 

7.50 

7.03 

3.92 

(3.71) 

D 

.5022 

(•5131) 

9.39 

(9.18) 

6.37 

3.98 

(3.76) 

.4 

.1988 

(.2094) 

7.70 

(7.59) 

4.09 

(3.95) 

.5 

.0909 

(.0980) 

.10 

6.45 

(6.40) 

4.4o 

3.78 

(3.69) 

I 


50 


TABLE  V  -  2 

a  »  .05/  «  26  -  26)  n^,  -  19 


\ 

OC 

ASN 

ELS 

Sl» 

A 

MPRT 

8PRT* 

IIPRT 

SPOT* 

14PFT 

0 

1. 

(1.) 

1. 

8. 

0. 

(0.) 

.1 

.9646 

(-9663) 

.95 

11.12 

3.72 

(3.53) 

.2 

.6975 

(.7059) 

15.08 

(14.60) 

5.21 

(5.07) 

s 

.5167 

(.5278) 

.50 

14.64 

(14.30) 

15.47 

11.85 

5.73 

(5.51) 

D 

.4711 

(.4826) 

14.43 

(Hf.l5) 

10.68 

3.84 

(5.6O) 

.4 

.1302 

(.1597) 

11.00 

(11.04) 

5.95 

(5.65) 

.5 

If 

.05 

8.72 

(8.89) 

6.57 

5.20 

(4.95) 
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TABIE  V  -  5 

a  -  .01,  »  47  -  4?)  •  38 


1. 


1. 

do 

.9995 

(.9993) 

.8084 

(.8151) 

.5066 

(.5175) 

.4247 

(.4358) 

.0185 

(.0205) 

.00076 

(.00092 


.999 


.50 


.001 
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ttB  V  -  4 

78  (0^-76)  n,.63 


Am 

HIS 

sm 

MPRT 

8PHT» 

MPRT 

27. 

(25.) 

18. 

0. 

(0.) 

56.42 

(34.96) 

28.83 

7.05 

(7.11) 

50.13 

(49.02) 

12.45 

(12.45) 

50.26 

(49.85) 

74.10 

45.24 

14.58 

(14.05) 

49.21 

(49.00) 

40.22 

15.20 

(14.57) 

29.98 

(30.84) 

13.90 

(13.60) 

21.31 

(22.14) 

17.03 

9.80 

(9.76) 
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A.k.2.  TABlfS  VI  -  1.  2.  k. 

\  *  1.  against  \  •  2. 
8  .  l.W^27  D  -  1.^570 


TABU!  VI  -  1 

a  -  .1,  Dq  -  17>  (o^  ■  U)»  Hp  - 


\ 

OC 

_ 

Asn 

HZ£ 

SEN 

A 

SPRT# 

IIPOT 

SPRT* 

MPRT 

0 

1. 

(1.) 

1. 

5. 

(3.) 

3. 

0. 

(0.) 

1. 

.9236 

(.9299) 

.90 

7.33 

(7.24) 

5.80 

2.66 

(2.85) 

1.35 

.6150 

(.6275) 

COCO 

?.46 

(3.44) 

• 

•5095 

(.5205) 

.50 

8.4l 

(8.45) 

6.97 

6.45 

3.55 

(5.51) 

D 

.4916 

(.5040) 

8.56 

(8.44) 

6.21 

3.56 

(5.52) 

1.75 

.2104 

(.2195) 

7.55 

(7.52) 

3.58 

(3.53) 

2. 

0617 

(.0867) 

.10 

6.16 

(6.56) 

4.61 

3.27 

(3.24) 
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Timi£  VI  -  2 

a  ■  .05,  -  26f  (a^  «  26),  1^  ■  I6 


\ 

OC 

ASN 

HIS 

SW 

MPRP 

SPRT* 

MPOT 

SPRT^^ 

MPRT 

0 

1. 

(1.) 

1. 

5. 

(5.) 

3. 

1. 

.9620 

(.9635) 

.95 

10.58 

(10.46) 

8.66 

1.35 

.6i»o6 

(.6if62) 

13.16 

(33.13) 

IBB 

s 

.5100 

(.5140) 

.50 

13.17 

(13.13) 

12.51 

10.83 

5.04 

(5.06) 

D 

.4879 

(.4936) 

33.09 

(33.10) 

10.42 

5.07 

(5.08) 

1.75 

.1573 

(.1607) 

11.07 

(11.15) 

5.12 

(5.08) 

2. 

.0413 

(.0427) 

.05 

8.80 

(8.91) 

6.88 

4.47 

(4.41) 
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TABLE  VI  -  3 

a  a  .01,  »  44,  ■  44),  aj  ■  32 


\ 

OC 

ASB 

HI£ 

SEN 

MPRT 

8PRT* 

MPRT 

SPRT* 

MPRT 

0 

1. 

(1.) 

1. 

8. 

(8.) 

5. 

0. 

(0.) 

1. 

.9922 

(.9926) 

.99 

-.J  cx> 

14.68 

5.67 

(5.70) 

1.35 

.6824 

(.6875) 

25.57 

(25.39) 

8.05 

(8.08) 

s 

.5066 

(.5083) 

.50 

25.85 

(25.60) 

30.46 

22.51 

8.30 

(8.34) 

D 

.4743 

(.4802) 

25.58 

(25.52) 

21.56 

8.43 

(8.39) 

1.75 

.0792 

(.0820) 

19.85 

(20.00) 

8.44 

(8.37) 

2, 

.0079 

(.0084) 

.01 

14.58 

(14.76) 

11.66 

6.50 

(6.50) 
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TABIE  VI  -  4 

a  -  .001,  -  71i  (Ogg-“"7l)>  nj  »  57 


\ 

OC 

ASR 

RIB 

sm 

msm 

8FRr« 

MPRT 

SPB®# 

MPRP 

0 

1. 

(1.) 

1. 

12. 

(I2O 

7. 

0. 

(0.) 

1. 

•9992 

(.99») 

.999 

28.44 

(27.90) 

22.46 

7.40 

(7.42) 

1.35 

.7315 

(.7358) 

44.73 

(44.41) 

12.07 

(12.15) 

8 

.5016 

(.5069) 

.50 

>5.1»3 

(45.33) 

68.82 

41.17 

12.58 

(12.52) 

D 

.4648 

(.4701) 

45.24 

(45.18) 

39.31 

32.71 

(12.63) 

1.75 

.0320 

(.C933) 

32.06 

(32.40) 

12.34 

(12.24) 

2. 

.00076 

(.00083) 

.001 

22.50 

(22.84) 

17.84 

ODCD 
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APPENDIX  V 

A. 5*  TABLE  VII  :  Characteristics  of  the  Poisson  liPRT 

»  D),  SPRTq#  SPRT  approximation 

and  FSST 


TABIES  VII  -  1,  2,  3>  ^ 

X.  >  .01  against  \  ■>  .1 

8  «  .03909  B  «  .01^299 


A.5.1-  TABIE  VII  -  1 

a  =  .1,  n^  -  58,  rij,  B  39,  n^  «  22.20,  n^^  »  19.08 


\ 

OC 

ASN 

SDN 

IIPHT 

SPHT^ 

0 

PSST 

IIPRT 

SPRT^ 

IIPHT 

SPRTq 

.0 

1. 

1. 

1. 

29. 

25. 

0. 

0. 

.0025 

.9958 

.9955 

1." 

29.98 

26.36 

3.84 

6.03 

.006 

.9776 

.9620 

.98 

31.04 

27.98 

5.60 

8.90 

.01 

.9432 

.9566 

.90* 

.94 

31. 8U 

29.44 

26.57* 

6.86 

10.89 

.015 

.8858 

.9116 

.88 

52.36 

30.72 

8.03 

12.53 

.02 

.8187 

.8553 

.81 

52.43 

51.45 

8.98 

13.68 

.025 

.7469 

.7916 

.75 

32.16 

31.73 

9.80 

14.52 

.03 

.6743 

.7240 

.67 

31.64 

31.63 

10.50 

15.17 

.035 

.6033 

.6555 

.60 

30.94 

31.24 

11.08 

15.66 

8 

.5479 

.6005 

.50* 

.55 

30.26 

30.74 

23.30* 

11.48 

15.56 

D 

.4976 

.5495 

.50 

29.36 

30.16 

11.79 

16.17 
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A.5«l-  TABLE  VII  -  1  (continued) 

a  »  .1,  -  58,  Up  -  39,  »  22.20,  -  19.O8 


A.5«2. 


TABLE  TEI  -  2 


a  ■  .05,  a  85,  Up  ■  63,  »  5r>53,  =  52.03 
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A.5.2.  TABLE  VII  -  2  (continued) 

a  *  .05,  n^  -  85,  nj,  ■  63,  •  37.55,  =  32.03 


X 

OC 

ASH 

SDH 

MPET 

SPRTq 

PSST 

IIPPT 

Hg 

IIPBT 

my 

.025 

.7703 

.8076 

.79 

48.85 

48.22 

16.13 

22.83 

.03 

.6844 

.7268 

.71 

48.23 

48.57 

17.54 

23.95 

.035 

.5978 

.6423 

.62 

47.08 

48.15 

18.74 

24.81 

s 

.5292 

.5738 

.55 

45.86 

47.54 

19.54 

25.32 

D 

.4672 

.50* 

.5106 

.49 

44.52 

41.84# 

46.26 

20.13 

25.67 

■  55 

.3054 

.3418 

.33 

39.82 

41.81 

21.12 

25.89 

.07 

.1681 

.1944 

.19 

35.86 

35.51 

20.80 

24.65 

.08 

.1097 

.1307 

.11 

30.27 

51.57 

19.97 

23.20 

.1 

.0448 

.0581 

.05 

24.35 

25.04 

17.68 

19.79 

.11 

.0282 

.05* 

.0388 

.03 

21.96 

18.95# 

22.46 

16.45 

18.09 

TABLE  VII  -  3 


a  B  .01,  n^  B  l43,  n^,  b  117,  b  78.39,  Hjj  ■  66.26 


X 

OC 

1  ASH 

L- _ _ _ 

SDH 

MPRT 

SPRT^ 

0 

SPBT^ 

0 

IIPRT 

SPRT 

.0 

1. 

1. 

l.„. 

67. 

52. 

0, 

0. 

.0025 

.9999 

.9999 

1.” 

71.69 

55.47 

6.59 

9.97 

.006 

.9991 

.9991 

1." 

75.75 

61.06 

10.96 

17.30 

.01 

.9931 

.9938 

.99 

80.62 

68.19 

14.95 

24.19 
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A. 5. 3*  Table  vii  -  3  (contin\^'>d) 

a  »  .01,  »  1^3,  ny  ■  II7,  -  78.39^  rxj^ 


oc 

ABS 

MPBT 

SPRTq 

FSST 

IIPRT  1 

SPRP^ 

0 

.01 

.99» 

67.25^^ 

.015 

.9681 

.9717 

.97 

86.39 

77.28 

.02 

.91^6 

.9223 

.91 

90.92 

1 

85.29 

.025 

.8315 

.8428 

.83 

93.55 

91.18 

.03 

.7264 

.7397 

.72 

94.05 

94.43 

.035 

.6107 

.6245 

.61 

92.59 

95.02 

.5164 

.5296 

.52 

90.19 

93.77 

.50» 

101.89* 

.4311 

.4435 

.44 

87.14 

91.39 

.055 

.2233 

.2327 

.25 

75.39 

79.66 

.07 

.0832 

1 

.0900 

.09 

60.69 

62.90 

.08 

.0401 

.0456 

.05 

52.52 

53.27 

.1 

.0064 

.0115 

.01 

40.46 

39.21 

.01* 

32.11* 

.U 

.0037 

.0059 

.05 

36.10 

34.31 

«  66.26 


SDIf 

MPRT 

SPBT 

0 

16.90 

30.68 

22.02 

34.97 

24.78 

37.73 

27.50 

39.75 

30.11 

41.49 

32.00 

42.78 

53.48 

43.82 

55.47 

44.99 

33.36 

41.61 

30.55 

37.56 

24.66 

28.99 

22.15 

25.31 

A.5-^. 


TABIE  VII  -  4 


61 


} 

1 


a  «  .001,  -  243,  Qy  »  212,  »  l43.97>  »  120.84 


\ 

oc 

SDIT 

MPRT 

PSOT 

MPRT 

SPRT^ 

0 

MPRT 

SPRTg 

.0 

1. 

1. 

1. 

110. 

77. 

0 

0 

.0025 

.9999 

.9999 

1" 

114.0 

82.33 

7.92 

12.50 

.006 

.9999 

.9999 

iT 

120.3 

91.07 

13.55 

22.54 

.01 

.9993 

.999^* 

iT 

128.5 

103.3 

19.41 

33.76 

.999* 

102. 9*^ 

.015 

.9917 

.9937 

.99 

140.0 

121.9 

26.13 

47.39 

.02 

.9610 

.9685 

.97 

151.3 

142.1 

31.53 

57.45 

.025 

.8903 

.9050 

.91 

160.2 

160.2 

35.75 

62.96 

.05 

.7763 

.7958 

.81 

164.8 

172.7 

39.95 

65.65 

.035 

.6331 

.6529 

.67 

164.0 

177.8 

44.78 

67.97 

s 

.5102 

.5272 

.55 

159.9 

176.7 

48.78 

69.61 

.50» 

230.2# 

D 

.3992 

.4124 

153.5 

171.1 

52.05 

73.07 

.055 

.1526 

.1564 

.18 

127.3 

140.9 

55.50 

76.19 

.07 

.0327 

.0338 

.04 

96.64 

100.8 

43.71 

65.68 

.08 

.0101 

.0110 

.01 

81. 4l 

81.01 

42.26 

55.22 

.1 

.00077 

.0012 

61.05 

56.20 

31.68 

37.96 

.001# 

49.15* 

.11 

.00020 

.00044 

48.51 

27.91 

32.07 
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APPEUDIX  VI 

A.  6.  TABLE  VIII;  Characteristics  of  the  Polason  IIPBT 

D)>  aPRT^i  SPOT  approxlmtlon  and 

ISOT 

TABUS  VIII  -  1,  2,  3, 

\  >=  2.  against  X  «  4. 
s  «  2.88^  D  «  2.914o 


A. 6.1.  TABLE  VIII  -  1 


a  ■  .1,  *^8  "  5'22>  n^j  »  3*10 


\ 

OC 

ASN 

sm 

A 

MPOT 

SPHT^ 

0 

PSST 

MPOT 

SPOT 

0 

MPOT 

SPOT^ 

0 

.0 

1. 

1. 

1. 

2. 

2. 

a 

0. 

1.5 

.9925 

.9952 

.99 

3.27 

3.08 

1.05 

1.34 

1.75 

.9758 

.9758 

.98 

3.62 

5.56 

1.27 

1.76 

2.0 

.9298 

.9410 

.94 

4.00 

4.12 

1.48 

2.19 

2.25 

.8U9I 

.90* 

.8658 

.86 

4.35 

2.90* 

4.68 

1.64 

2.54 

2.5 

.7507 

.7462 

.75 

4.56 

5.11 

1.75 

2.78 

2.65 

.6463 

.6572 

.67 

4.63 

5.27 

1.80 

2.87 

2.75 

.5871 

.5958 

.61 

4.64 

5.32 

1.83 

2.91 

s 

.5062 

.5065 

.54 

4.62 

5.32 

1.86 

2.95 

D 

.4893 

.50* 

.4883 

.52 

4.61 

3.48<f 

5.51 

1.86 

2.95 

5.25 

.5070 

.2938 

.54 

4.37 

4.98 

1.89 

2.93 

64 


A.6.1.  TABIE  VIII  -  1  (continued) 


a  «  .1,  n^  B  10,  n^B  6,  n^  «»  5-22,  n^^  *  J.IO 


\ 

OC 

ASIT 

sm 

UPRT 

aPBT^ 

0 

psar 

MPBT 

SPBT^ 

0 

MPBT 

SPBT^ 

1  0 

5-5 

.2016 

.1858 

.23 

4.08 

1.86 

2.80 

3.75 

.1254 

.1116 

.14 

3.77 

1.79 

2.60 

4. 

.0746 

.0649 

.09 

3.45 

3.64 

1.69 

2.34 

4.5 

.0240 

.1* 

.0211 

.03 

2.90 

2.30* 

2.90 

1.45 

1.82 

A.6.2.  TABIE  VIII  -  2 


a  ■  .05,  n^  -  13,  nj  «  8,  ng  -  -  5.21 


\ 

OC 

ASIT 

SEN 

MPRT 

SPBTq 

PSST* 

MPBT 

SPBT^ 

MPBT 

SPBTq 

.0 

1. 

1. 

1. 

3. 

2. 

0. 

0. 

1.5 

.9984 

.9984 

.99 

4.35 

3.76 

1.38 

1.65 

1.75 

.9912 

.9925 

.99 

4.96 

4.48 

1.72 

2.28 

2.0 

.9647 

.97U 

.95 

5.65 

5.42 

2.04 

2.95 

2.25 

.8973 

.95* 

.9118 

.86 

6.33 

4.33* 

6.45 

2.29 

3.50 

2.5 

.7740 

.7931^ 

.72 

6.85 

7.34 

2.46 

3.82 

2.65 

.6761 

.6944 

.62 

7.02 

7.68 

i 

2.53 

3.91 

2.75 

.6047 

.6206 

.55 

7.06 

1 

7.80 

2.58 

3.95 

8 

.5052 

.5166 

.50* 

.46 

7.03 

7.83 

6.25* 

2.63 

3.99 
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A.6.2.  table  VIII  -  2 


a  «  .05,  ■  13>  n*  «  8,  ■  5*^»  nj,  «  5*21 


\ 

oc 

ASIT 

SDIY 

A  ' 

MPRP 

SPBT^ 

0 

PSST* 

MPRT 

1 

SPRTq 

MPBT 

SPBT^ 

0 

D 

.i»8U3 

.4947 

.44 

7.01 

7.82 

2.65 

3.99 

3.25 

.2648 

.2644 

.25 

6.51 

7.20 

2.73 

3.9k 

3.5 

.1496 

.1466 

.15 

5.93 

6.42 

2.68 

3.75 

3.75 

.0772 

.0752 

.08 

5.30 

5.59 

2.54 

3.42 

i^. 

.0370 

.0369 

.05 

4.72 

4.83 

2.33 

3.00 

.05* 

3.44* 

4.5 

.0073 

.0088 

.01 

3.78 

3.70 

1.87 

2.20 

A.  6.3.  TABIBVIII-3 

a  ■  .01,  ■  22,  Hy  ■  16,  ■  11.25»  Op  ■  10. 78 


oc 

ASH 

sw 

m 

inQ^^iiiiiiiii 

SPRT^ 

PSST# 

MPRT 

SPRT^j 

MPRT 

SPRPq 

.0 

1. 

1. 

1.. 

4. 

3. 

0. 

0. 

1.5 

QOOO 

.9999 

!,• 

7.05 

5.61 

1.67 

2.11 

1.75 

.9993 

.9995 

1." 

8.06 

6.79 

2.23 

2.96 

2. 

.9928 

.99*H 

.99 

9.39 

8.51 

2.91 

4.20 

2.25 

.9567 

.99* 

.9641 

.94 

11.00 

7.33» 

10.84 

3.57 

5.49 

2.5 

.8442 

.8601 

.79 

12.54 

13.31 

3.97 

6.27 

A.6.5. 


TABIiB  VIII  -  3 


a  *  .01,  =  22,  Hj,  B  l6,  « 


\ 

oc 

AS^ 

IIPRT 

SPBTq 

FSST# 

MPRP 

2.65 

.7285 

.7^58 

.66 

13.16 

2.75 

.6356 

.6511 

.57 

13.37 

s 

.5008 

.5112 

.50* 

.55 

13.35 

D 

.4723 

.4814 

.42 

33.30 

3.25 

.1908 

.1880 

.17 

11.88 

3.5 

.0757 

.0723 

.07 

10.32 

3.75 

.0250 

.0237 

.03 

8.84 

4. 

.0071 

.0072 

.01# 

.003 

7.62 

4.5 

.00042 

.00075 

.001 

5.93 

A. 6. 4.  TABLE  VIII  -  4 


a  »  .001,  ■  36,  iiy  ■  28, 


V 

oc 

m 

SPBTq 

PSST# 

MPBT 

— 

1. 

1. 

1. 

6. 

.9999 

.9999 

r 

10.83 

.9999 

.9999 

r 

12.42 

2. 

.9993 

.9994 

.999* 

r 

14.55 

66 
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(continued) 

11.25,  rjjj  =  10.78 


r 

sm 

SPRP^ 

0 

MPRT 

14.41 

4.10 

6.43 

14.84 

4.17 

6.49 

14.98 

15.23# 

4.30 

6.56 

14.95 

4.32 

6.58 

13.11 

4.52 

6.65 

10.93 

4.31 

6.18 

8.89 

3.84 

5.30 

7.29 

5.83# 

3.29 

4.32 

_J 

5.27 

2.39 

2.84 

m 

20.59,  Hjj  - 

19.65 

1 

SDIf 

SPBT^ 

MPRT 

SPRT^ 

4; 

0 

0 

7.86 

2.10 

2.49 

9.58 

2.79 

3.65 

12.25 

11.23# 

3.75 

5.53 

A. 6. 4. 


TABIE  VIII  -  4  (continued) 


a  B  .001#  n^  •  36#  Up  ■  28#  n^  ■  20. 59^  Djj  i 


/\ 

MPHT 

SPRT^ 

.  0  ! 

1 

PSST* 

MPHT 

SPRT 

0 

2.25 

.9871 

.9898 

.99 

17.1^5 

16.51 

2.5 

.9060 

.9151 

.91 

20.85 

22.15 

2.65 

.7856 

.79^^3 

.79 

22.51 

25.15 

2.75 

.6738 

.6783 

.68 

23.16 

26.46 

s 

.5005 

.4963 

23.26 

26.99 

34.41* 

D 

.4634 

.457!^ 

.44 

23.16 

26.91 

3.25 

.1269 

.1143 

.14 

19.79 

22.17 

3.5 

.0307 

.0255 

.04 

16.43 

17.10 

3.75 

.0053 

.0044 

.01 

15.64 

15.09 

4. 

.00069 

.OOOTfk 

.001 

11.57 

10.38 

.001* 

8.92* 

'*.5 

.00071 

1 

.00034 

.ooos. 

5  8.89 

7.32 

67 


«  19.65 


OT 


MPRT 

SPRT^ 

4.97 

8.13 

5.89 

10.01 

6.08 

10.22 

6.16 

10.15 

6.38 

10.18 

6.44 

10.21 

6.88 

10.64 

6.24 

9.i^l 

5.19 

7.37 

4.24 

5.57 

2.99 


3.46 


A. 6.5.  ASN  CURVCS 


■I  ■  001 


Figurt  3.  Poimon  ASN 
H| :  A  -2.  against  H2:  A  -  4. 

.  5^2.8853901  .A-2.SI39S9 
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COMPILER  PROGRAM 
A, 7,  K.  Fukushima.  MPRT-A 


Y  2400 


^thli 
the 


kfbr 
keaoa 


Luea 


B957-ftheM 
n998«fare  thi 
n959-^  uppe| 
n9 

a^l-siomr 
n942-iboiaidai 
n943«^M^ 
n944>4theae  $ 
n945-iarc  thi 
n946«Ae  aceei 
n947-lptaee^ 
n94^  probef 
H9494bllltli 
n9$0>$e8  forf 


Z  2400 
f 
f 
t 
t 
t 
t 
t 
t 
t 
t 
t 
f 
t 
t 
f 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 

r 

f 

t 

t 

t 

t 

t 

t 

t 

t 

t 


8  0190  V  0200 


I 


COMPILER  IROCRAM  (cont.) 


a951*$  eaoh 
n958-f trial 
a9?3*?tte8e 
a95M 

a955-^  reje^ 
a9?o*ietloa 
b957- 

a^^llltle^ 

a96c . 
b96i< 
a962«4thl8  li 

B9o4«^proD«  9 
a965*$of  aeei 
a966>|eptane4 
B967-$e  I 
a96^thla  li 
a969^B  the  i 
B970-fprob.  I 


a903- 
a964^nwiberi 
a^«i  of  eai 
B^xiiivlee  f 
a9674tble  li 
the  i 


alOOO-ie  evani 
BlOCn-ilatlTei 
al002-i  poleei 


f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 
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COtPILER  FROCStAM  (cont.) 

ii1003>^  proA 
nlOO^ib.  ftrol 
nnOO!$-n  seroi 
n1006«4the8e  % 
n1 007-$are 
ii100e>d 

nl 009*9latlve^ 

Bl(no«i  polsi 
n1(ni«^ii  pi 
ol012>ib.  yjcp 
n1 01 34^0  lnf$ 
n10l4>ilnlty  $ 
n1015-$thls  li 
n10l6«$8  the  4 
ofi 

a10l8-i  trlali 

n1 019*4  pr*et4 

nl  020>>4ically4 
n1021>4^hl8  l4 
nl022>48  the  4 

n1023*4**^*  ^ 

nl02i^f  numbi 
nl023sMr  of  i 
nl  026«4mbip1o4 

nl 027*4*  4 


71 


f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

t 


0001  yl,...,y4  y^,.,.,y36  fukoehim  iapnt 

011  6  tyl  ty2  ty3  ty4  f 
013®  ty^  ty29  ty»  f 
0117  ty31  ty32  ty33  tyjk  t 

0118  ty35  ty36  ty37  ty38  f 

*l«(*0l8,y3*  -  *01e*yl*)  f 

*2»(*Ol8,y3*  -  *0l8,y8»)  f 

*3-j*0l8,y2*  -  *0l8,y1*)  f 

*4*(*0l8,2xy4*)  f 

y5>sU/B2  f 
y6-(y3-y2)/«2  f 

y7.((-i)»4y.3  f 

y0-(y2-yl  )/s3  f 

0080  ty5  ty6  ty7  ty6  f 

0079  ntn903  ntn913  f 

n90e-y7-y949  f 

tn902  f 

atn9l4  nta9l9  f 

y9-(y7-y5)/(y6-ye)  f 
0078  ty9  f 

0077  nta920  atn923  t 

ylO-y54y63cy9  f 

y11-y7-1y6jty9  f 

0076  tyio  tyii  f 

0073  atn924  ata931  f 

0074  atn932  atn936  f 

nl-O  f 

0002  y(1000Hil)-y5^y6xnl  f 
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COMPILER  HiOORAM  (cont.) 


Ml  •yi+y&xnl  f 

g3  if  y(l  OOOfni )  w  0  f 

yOOOOhil )-(-!.)  f 

0003  njl Ohil  )-y(l 0004ti1 )  f 

a(30O4iil  )>zii1  f 

gU  If  zn1.n(500tn1)  u  0  f 

n(500fnl)-n(500tel)+l  f 

0004  yO-yO  f 

tn(50CMn1)  tal  tB(IOhil)  f 

m-nl+1  f 


g2  If  a(500«ti1.1).B(10ln1.1)  y  i  f 
ii900»nl-'|  f 

atn937  atn943  f 

0078  tn900  f 

0071  atn10l3  ata1080  f  The  end  of  Pert  A. 

y39-y2  t 

n901-0  f 

y2-y26  f 

g40  If  yp  u  0  f 

^  f 

0040  n9Gfl-1  f 

y2«y29  f 

g4l  If  yp  u  0  f 

85  f 

0041  B901-P  f 

y2-y30  f 

iflSO  If  yP  V  0  f 

85  f 

0180  B901-3  f 

y8-y5i  f 

0181  If  y8  u  0  f 
85  f 

0181  b901>4  f 

y8-y38  f 

88  If  yP  tt  0  f 

65  f 

0188  11901-5  f 

ya-y35  f 

tf23  if  yP  U  0  f 
85  f 

0183  b901-6  f 

y8-y34  f 

1^84  If  yp  u  0  f 
85  f 

0184  n901-7  f 

y8-y35  f 

0125  if  yP  u  0  f 
85  f 

0185  0901-8  f 

y2-y36  f 

8126  If  yP  u  0  f 
85  f 


73 


COMPILER  PROGRAM  (cont.) 

0126  n901-9  f 

y2-y37  f 

efI27  If  y2  u  0  f 

g5  f 

0127  n901-10  f 
y2«y30  f 
g42  if  y2  u  0  f 

85  f 

0005  y1000-(*028,(-l)xy2*)  f 

tyiooo  f 

7,n2,1,l,n902,  f 

0006  y12>n2  f 

y(1000«n2)-(y(1000ki2-1)xy2)/y12  f 

n2«y12  f 
ty(1000hi2)  ta2  f 

0007  yO»yO  f 

atn994  atn997  f  end  of  Part  B. 

*1000-yi0O0  f 

tziOOO  t(0.)  f 

9,n2,l  ,1 ,n902,  f 

0000  zOOOOHi2)*a(1000ta2-1)4y(1000«n2)  f 

t2(1000fn2)  ta2  f 

0009  yO=«yO  f 

atn998  atn1005  f  The  end  of  Part  C. 

*0«1  f 

t(1.)  t(0.)  f 

51 ,n2,l ,1 ,n902,  f 

0050  cii2>1.c(1000«ti2-1)  f 

tznP  tii2  f 
0051  yO«yO  f 

atn1006  atn10l4  f  The  end  of  Part  D« 

m-1  f 
m050=n501-1  t 

0011  12,n2,nl  1,1 ,111050,  f 

0012  y(100fn2)-y(l  000*02)  f 

8^3  If  mi  u  (-1)  f 

y2001-y(  100*011)  f 

ty2001  t(l.)  f 

f 

0015  y2001-0  f 

ty20Cn  t(1.)  f 

0014  12001-10501  f 

t(0.)  t(l.)  ti20Q1  f  The  end  of  Part  E. 
0015  nl-nl+1  f  Part  F 

0017  gZl  lfn(10tnl)v0  f  part  G 

y(2000*ol)»0  f 

t(0.)  tnl  f 

m051-n(500*nl)>l  f 

n10^-n(500*ol-l)  f 

n1 095-n(500*o1  -1  )-l  f 

20,n2,0,l  ,01095,  f 


0018 

0019 

00|6 

0020 


0091 

0092 

0093 

0094 

OOOfI 


0022 


0023 

6023 

0026 

0027 

0096 

0003 

0006 

0007 

0090 

0020 

0029 


COMPILER  PROCBIAM  (cont.) 


y(300h>2)-0  f 

lOfii3,0,1,n2,  t 

y( 300tii2)-y( 300lti2)4y(l  00«n3)xy(1  OOOMi2-a3) 
ty(3001fi2)  tal  ta2  f 

yO^^  f 

g26  If  ii(300hH  -1 )  II  a(900M1 )  f 

94,ii2«ii1  094,1  ,n1 031,  f 

y(300fii2)-0  f 

93,113,0,1  ,n1 095,  f 

y(300hi2)-y(3eo«fi2)-ly(l  00«tii3)xy(1  OCOIat  ■3) 
ty(300tfi2)  tal  ta2  f 

yO^i^  f 


82®  f  The  end  of  Part  H. 

g22  If  n(10Hi1)-B(10»fi1-1)  V  0  f 

y(200dhi1  )-0  f 

t(0.)  tal  f 

Bl034-n(10fti1)+1  f 

a1033>m(300Hi1)-1  f 

625  f 

y(2000kt1  )-0  f 

n1032«anotn1 -1)4-1  f 

Bl03>>«n0hi1)  f 

n1034>nn(Hn1)4-1  f 

n1035-B(500W1)-1  f 

23,a4,n1 032,1 ,a1 033,  f 

y(20004a1  )-y(2000hi1  )4y(lOO«d4)xB(1000«e(10M)HiA) 

ty(2000ki1)  tal  f 

ata944  ata952  f 

a1031-a(300«n1)-1  f 

Bl094«n(300*n1-1)  f 

Bl095«n(300fn1-1)-1  f 

o1032-n(10Hi1-1  )4-1  f 

90,n2,n1 034,1  ,Bl  093,  f 

y(300Hi2)-0  f 


27,B3,m  032,1  ,b2,  f 

y( 300t«i2)-y( 300hi2)4r(1 00*tt3)^1 000tflft-a5) 
M300tii2}  tal 035  ta2  ta3  f 
yO^  f 

if  B(300ta1-1)  u  a(300tn1)  f 

90,B2,a1 0^,1  ,a1 031,  f 

^300Hi2)-0  f 

07,b3,b1  052,1  ,a1 095,  f 

y(  300Hi2)-y( 300hi2)4r(1 00hi3)xy(1 000«b8-b3) 
ty(300Hi2)  tal  ta8  f 

S2?500kif)  f 

a1057-B(10hi1 -1)4-1  f 
Bl058-B(500tBl-l)-1  f 

y15-0  f 


29,b5,b1  057,1  ,Bl  050,  f 

yl5-yl54yO«>hi5)xi(B2-B5)  f 
*(2000tn1)»y15  f 
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f 


f 


f 


f 


f 
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COUPILES  FRO(»AM  (cont.) 
tnl  tz(2000»n1)  r 

atn9!53  ata96l  f  The  end  of  Fart  J. 

mo6i-nno4fn)+i  t 

n1062«n(30Ofn1)*1  f 

31  ,n2,n1 061,1, n106e«  f 

0031  y(l004i\2)-y(300«ti2)  f 

g32  If  n1  w  n900  f  Part  K 

fifl5  t 

0032  y20-0  f 

y21-0  f 

y22»0  f 

y23-0  f 

3U,n1,1,1,n900,  f 

0033  y20-y20«yf 2000«tt1 )  f 

y21-y2l+*(2000lin)  f 

y22-y22fn1x(y(2000tailHtf2000MI))  f 
0034  y23-y23hi1ai1x(y(2000«Bl)4a(200ehA))  f 
0103  ty20  tn900  ty2  f 

mo6  ata962  ata967  f 

0107  ty21  tn900  ty2  tyU  f 

0106  atn968  ata972  f 

0109  ty22  ty2  ty4  f 

OHIO  atn973  ata979  t 

y24-y23-(y22xy22)  f 

0111  ty24  ty2  ty4  f 

0112  atn960  ati^  f 

y25-(*06i,y24*)  f 

0113  ty25  f 

OIIH  ata1021  ata1027  f 

0115  atn967  atn993  f  The  end  of  Part  L. 

g4o  If  a901  u  0  f 

g4l  If  n90l  u  1  f 

ffl20  If  n901  u  2  f 

e^21  If  n901  u  3  t 

e(l22  If  n90l  u  4  f 

g\23  If  a901  u  5  f 

C(I24  If  n901  u  6  f 

cn25  If  n90l  u  7  f 

e^26  If  a901  u  8  f 

C(127  If  11901  u  9  f 

g42  If  n90l  u  10  f 

0042  y2«y39  f 

0043  21  tt 


The  end  of  Par  M. 
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